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The magneto-photoluminescence in modulation doped core-multishell nanowires is predicted as 
a function of photo-excitation intensity in non-perturbative transverse magnetic fields. We use a 
self-consistent field approach within the effective mass approximation to determine the photoexcited 
electron and hole populations, including the complex composition and anisotropic geometry of the 
nano-material. The evolution of the photoluminescence is analyzed as a function of i) photo¬ 
excitation power, a) magnetic field intensity. Hi) type of doping, and iv) anisotropy with respect to 
field orientation. 


I. INTRODUCTION 

In the last decade, semiconductor nanowires (NWs) 
have received increasing attention owing to their po¬ 
tential applications in nanotechnology.^ They represent 
promising candidates for a wide range of novel ultrafast 
electronic and optoelectronic nanodevices including high 
electron mobility transistors,^ photovoltaic cells,^ single¬ 
photon emitters,^ and lasers.^ Gritical steps have been 
taken in these directions in the last years. Long single¬ 
crystal, defect-free cores,selective radial doping,^ lat¬ 
eral overgrowth with high-quality interfaces^ have been 
successfully realized. Most importantly, complex radial 
modulation doped heterostructures are now being engi¬ 
neered in core-multishell More complex than 

their planar counterparts, radial heterostructures in these 
systems may host high-mobility electron/hole gases with 
an inhomogeneous localization in the section of the NW, 
which strongly depends on doping density and gate po¬ 
tentials. 

Magnetic states of such axial electron gases have been 
investigated in cylindrical systems^^^^^ or in wrapped 
heterojunctions with prismatic symmetry.In the 
Quantum Hall regime, the magnetic field competes 
with the self-consistent confinement field. Recent self- 
consistent field calculations of realistic GaAs/AlGaAs 
core-shell NWs in transverse fields^^ showed that at 
high carrier density the carrier-carrier repulsion results 
in complex magnetic bands, corresponding to inhomo¬ 
geneous charge localization at the prismatic hetero¬ 
interfaces of the nanomaterial. 

Intraband optical spectroscopy is a key component in 
the investigation of properties of single semiconductor 
NWs and NW heterostructures.In particular, 
photoluminescence (PL) has been successfully used in 
bare and core-multishell (GMS) NWs in order to analyze 
the dynamics of photo-excited electron-hole plasmas^^”^^ 


and the quantum confinement effects.Magneto-PL 
in GaAs/AlAs GMS NWs has been applied to undoped^^ 
and doped samples,where the magnetic field depen¬ 
dence of peaks ascribed to confined carrier emission has 
been examined. 

In this work, we consider magneto-PL in modula¬ 
tion doped GaAs/AlAs CMS NWs with an inversion or 
an accumulation layer at the inner GaAs/AlAs hetero¬ 
interface in a transverse magnetic field. Using a self- 
consistent envelope-function approach, we predict the 
magneto-PL at different excitation powers and identify 
non-linear effects due to filling and restructuring of the 
free-carrier gas. Furthermore, we investigate the effects 
of the type of doping and of the structure anisotropy on 
the magneto-PL. 

The article is organized as follows. In Sec. II, we de¬ 
scribe the physical and numerical modeling of the mod¬ 
ulation doped CMS NW used to evaluate the hole and 
electron states and from these the PL spectra. The re¬ 
sults obtained from our numerical calculations are dis¬ 
cussed in Sec. III. First, the prototype nanostructure 
under investigation and its relevant physical parameters 
are presented in subsection III A. Then, subsection III B 
illustrates predictions for the p- and n-doped system at 
various regimes of photo-excitation power. Anisotropy 
with respect to field orientation is examined in subsec¬ 
tion IIIC. We summarize our results in Sec. IV. 


II. THEORETICAL APPROACH 

Electron and hole states are determined by a self- 
consistent Schrodinger-Poisson approach.Within 
a single parabolic band effective-mass approximation, the 
Hamiltonians for holes {h, upper sign) and electrons (e, 
lower sign) in an external magnetic field^^ are given by 
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respectively, where R=(x,^,z) represents a 3D coordi¬ 
nate, P is the conjugate momentum operator, q is the 
positive unit charge, A(R) is the magnetic vector poten¬ 
tial, (R) is the material-dependent effective mass of 
electrons (holes), and ^y(R), EcCR), and ^(R) are the 
local valence-band edge, the local conduction-band edge, 
and the electrostatic potential generated by the free car¬ 
riers and the dopants, respectively. Electron-hole interac¬ 
tions are taken into account within a mean-field approx¬ 
imation through the self-consistent potential E(R), i.e., 
’excitonic’ effects are not included in our calculations, as 
they are expected to be negligible at densities larger than 
10^ electron/cm^^’^^^^ 



FIG. 1. Sketch of the GaAs/AlAs core-multi shell NW in¬ 
vestigated in this work. The GaAs core is overgrown with 
AlAs/GaAs/AlAs layers forming a narrow quantum well and 
a GaAs capping layer. The GaAs quantum well serves as 
a doping layer with a homogeneously distributed density of 
donors or acceptors ud or tia), respectively. The transverse 
magnetic field forms an angle 0 with the x-axis. 


The transverse magnetic field is described in the Lan¬ 
dau gauge A(R) = B{y cos 0 — x sin 0)z^ where 0 denotes 
the angle between the field and the x-axis (see Fig. 1). 
We shall consider two possible orientations of the field: 
i) perpendicular to the facets {0 = 7r/2) and ii) along the 
maximal diameter {0 = tt/S). 


Since the field does not break the translational in¬ 
variance, the 3D wavefunctions can be factorized as 
^n^^R) = where r = {x,y) and k the 

wavevector along the nanowire axis. Then the 2D wave- 
functions (r) are solutions of the Schrodinger equa¬ 
tions 
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where uJc,e{h){'^) = ||)(^) (electron) 

cyclotron frequency, Ib = ^/hl{qW) is the magnetic 

length, and is the energy of the hole (electron) 

state. Note that distinguish between the hole effective 
mass in the section of the CSNW, and the hole 

effective mass along the wire, || ’ tske into consid¬ 
eration the strong mass anisotropy of GaAs. The electron 
effective mass, m*, is assumed to be isotropic. 

Similarly to the familiar case of a planar Hall bar, in 
the above equations the effect of the magnetic field is 
represented by an harmonic potential displaced by kl^ 
from the nanowire axis. Therefore, in general, is 

not parabolic in k at large magnetic fields. Note that 
Eqs. 2) are coupled by the electrostatic field V (r) which is 
determined by both electron and hole densities, as shown 
below. 


Equations (2) are numerically solved by the box inte¬ 
gration method on a triangular grid with hexagonal el¬ 
ements. At each magnetic fields the equations are 
solved on a grid of wavevectors in [0,/cmax]-^^ The k- 
domain is discretized in steps Ak = 0.075 nm“^, while 
kmax is taken sufficiently larger than the Eermi wave vec¬ 
tor that its occupation is negligible. The hole and elec¬ 
tron densities are evaluated at temperature T from 
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where the factor 2 accounts for spin degeneracy, and 
f{x) = 1/(1 -h e^) is the Eermi function. Ep indicates 
the Eermi level which is determined by surface states 
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and, therefore, by the type of doping. The charge densi¬ 
ties in Eqs. (3) depend on the fully-ionized dopants, i.e., 
photo-excited carriers are not included here. As shown 
elsewhere,the symmetry and localization of and 
Ug are strongly influenced by the level of doping, due to 
competing energy scales.The linear free charge den¬ 
sity per unit length provided by the dopants can 

finally be calculated from 

<e(ft) = [ K(h)i^’EF)dr. (4) 

J A 

Its dependence from Ep omitted for brevity. 

The laser field also excites electron-hole pairs with a 
linear density If the photo-excited density is not neg¬ 

ligible with respect to the free charge induced by dopants, 
the self-consistent field should take these charges into ac¬ 
count. The total electron (hole) linear density is 

n\,e(h) =nl^ih)+A'"- ( 5 ) 

In order take into account the photoexcited charge, we 
shift the Fermi level Ep to E^^^ (in general Ep ^ Ep) 
until 
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where the product 

of the two Fermi distributions describes the occupation 
of the states involved in the photoemission, and the last 
term accounts for a Lorentzian broadening of the peaks 
due to incoherent processes by means of a phenomeno¬ 
logical parameter 7 . As usual for interband transitions 
at optical frequencies, only vertical transitions are taken 
into account. 


This amounts to the charge neutrality condition in an 
undoped sample. In practice, we obtain E^p^^ from (3)- 
( 6 ) by means of the bisection method. 

Once the total electron and hole densities are com¬ 
puted, the electrostatic potential V{t) is obtained by 
solving the Poisson equation 
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where e^(r) and eo describe the position-dependent di¬ 
electric constant and the vacuum permittivity, while 
^D(A)(i*) denotes the ionized donor (acceptor) density. 
Dirichlet boundary conditions are used to solve Eq. (7) 
with the potential on the domain boundaries set to zero. 

The electrostatic potential is inserted into Eqs. ( 2 ) and 
the whole procedure is iterated until self-consistency is 
achieved, that is, the relative variation of both the elec¬ 
tron and hole density is < 0.01 at each point of the sim¬ 
ulation domain. 

The PL spectra is computed under the assumption that 
carriers are thermalized in the lowest electronic levels by 
fast non-radiative processes. Thus, the intensity of PL 
signal is given by^^ 
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TABLE 1. Material parameters used in the numerical calcu¬ 
lations.^^ The relative effective masses in the axial {z) and 
transverse (r) directions are given in units of the free electron 
mass. 


III. RESULTS 


A. Structure details 

We simulate a prototype NW with the same composi¬ 
tional parameters as in Ref. 12 . A GaAs core of diameter 
60 nm (facet to facet) is surrounded by a 3 nm/1 nm/3 
nm AlAs/GaAs/AlAs multilayer and a 13 nm GaAs cap¬ 
ping layer. Donors or acceptors are uniformly distributed 
in the narrow GaAs quantum well. For p-doping we as¬ 
sume a Fermi level pining at 400 meV above the valence 
band edge,^^ while for n-doping we assume mid gap pin¬ 
ning. The PL phenomenological linewidth 7 = 0.5 meV 
and all calculations have been performed at temperature 


T = 1.8 K. 

Material parameters are reported in Tab. 1. Hole 
masses are chosen as those of a planar quantum well 
grown along the [ 110 ] direction, which corresponds to the 
radial quantum well of the GSNW structure. Appropri¬ 
ate expressions of the effective masses along the quanti¬ 
zation axis [ 110 ] of the quantum well and the in-plane di¬ 
rection (here corresponding to the in-wire direction [ 111 ]) 
are obtained, e.g., in Ref. 36 from the diagonalization of 
the Luttinger Hamiltonian. Note that the ’heavy’ mass 
along the [ 110 ] direction is taken as isotropic in the sec¬ 
tion of the GSNW, why a much lighter mass describes 
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the heavy hole dispersion along the wire. 

Note that in our procedure we do not include light hole 
states. We have checked that light-hole subbands (eval¬ 
uated within a parabolic approximation with the light- 
hole mass in the [110] growth direction^^) lie well be¬ 
low than the Fermi level. Clearly, light-hole states might 
still contribute to heavy-hole states via spin-orbit cou¬ 
pling which is neglected in the parabolic approximation 
assumed in this work. Unfortunately, spin-orbit coupled 
states have not been investigated for realistic models of 
CSNWs. However, from previous calculations in semi¬ 
conductor quantum wires, we may expect these effects to 
be small for these weakly confined states.^^ 

Note that the narrow GaAs QW acts here as a dop¬ 
ing layer. Due to its narrow thickness, no quantum state 
is hosted by the QW itself, and all free charge accumu¬ 
lates in an inversion or accumulation layer, depending on 
the type of doping, in the core or at the inner hetero¬ 


interface, as we shall discuss below. 

B. Fan diagrams 

We examine the evolution of the magneto-PL with a 
transverse magnetic field for p- and n-doped samples. We 
shall consider three levels of the photo-excited charge, 
namely, negligible, lower or comparable to the free-charge 
density induced by a typical doping level. 

1. p-doping 

We first consider a p-doped sample with ua = 2.5 x 
10^^cm“^, which corresponds to a linear free-hole density 
Kh = 8-5 X lO^cm We first consider a low excitation 

density = 10^ cm“^ <C 




-60 -30 0 30 60 -40 -20 0 20 40 -30 -15 0 15 30 

k Ig(nin) k Ig(nm) k Ig(nm) 


Intensity(arb. imits) 

0 0.25 0.5 0.75 1 

1560 

1550 
> 1540 

s 

PL) 

1530 

1520 

0 5 10 15 20 

B(T) 



FIG. 2. Magnetic states and magneto-PL calculated for ua = 2.5 x 10^^ cm“^ and = 10^ cm“^. Here B is normal to the 
facets (0 = 7r/2). Left: electron (top) and hole (bottom) subbands vs kl% at selected fields B as indicated. The horizontal lines 
indicate the energy of the electron/hole Fermi level. The hexagonal insets show the electron (top) and hole (bottom) charge 
density at H = 2 T. Right: calculated PL spectra. The inset displays the squared envelope functions for the electron (left) and 
hole (right) state whose recombination yields the strong peak at F = 1523 meV when F = 0. The regular intensity oscillations 
are a graphic artifact due to the course magnetic field grid. 


In Fig. 2, left panel, we show the dispersion of the 
magnetic levels for holes and electrons evaluated at se¬ 
lected B values. The evolution of the magnetic subbands 
for conduction electrons resembles that of Landau levels 
and edge states in planar Hall bars.^^ Indeed, as shown 
in the inset, the electron charge is localized in the core, 
since is very small and it is not able to bend the 
conduction edge. Therefore, electron states are similar 
to those of a QW in a vertical field, here represented by 
the NW core, and bend up at increasing k, which corre¬ 


spond to edge states confined by the field on the lateral 
boundaries of the GaAs core.^^ 

This is at difference with the dispersion of holes, which 
becomes flat at lower fields and develops non monotonous 
behavior at large k. Such a behavior is not due to the 
strong anisotropy of the hole mass but to the emergence 
of peculiar localized states.Indeed, holes, which here 
are mainly generated by doping so that riph ^ cre¬ 
ate an accumulation layer near the hetero-interface, with 
an approximately tubular geometry, as shown in the in- 
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set. As discussed in Ref. 18, at sufficiently large fields 
the energy dispersion curves develop minima which corre¬ 
spond to states localized at the flanks of the tube, where 
the vertical component of the field with respect to the 
charge layer (and the associated confinement energy) is 
vanishing. Note from Fig. 2(left) that several hole lev¬ 
els are occupied, while electrons only occupy the lowest 
magnetic subband. The number of occupied subbands 
in the hole states decrease with increasing field. At in¬ 
termediate fields {B = 10 T), only 4 hole subbands are 
occupied and at 22 T only two hole magnetic subbands 
he below or close to the hole Fermi level. 

The evolution of the ground state n = 1 and k = 0 for 
electrons and holes is shown in Fig. 3 at selected magnetic 
field intensities. While the electron state is delocalized 
over the core, the hole state is localized at the top and 
bottom of the structure (here the field is in the vertical 
direction). Clearly, the states extend laterally by about 
Ib- 

The evolution of the calculated PL spectra with field is 
shown in the fan diagrams of Fig. 2, right panel. At zero 
field, the spectrum exhibits two peaks. The lowest in en¬ 
ergy, close to the band-gap of bulk GaAs at 1519 meV, is 
due to the electron and hole (e-h) ground state, coupling 
the states shown in Fig. 3. Due to the small overlap, this 
peak is weak. The second, much stronger peak originates 
from the recombination of the electron ground state and 
the hole lowest excited state spreading over the center of 
the GaAs core, hence the large oscillator strength (see 
the inset in the right panel of Fig. 2). 

In the PL spectra in Fig. 2 we can distinguish two 
different field regimes. At low fields, the peaks undergo 
a quadratic diamagnetic shift and progressively reduce 
their intensity. The stronger peak disappears at 5 ^ 10 
T, due to the decreasing occupation of the hole excited 
state (n=3). Indeed, even though the global hole density 
is almost independent from 5, as the field increases the 
subbands flatten and accommodate a larger number of 
carriers in the lowest subband with k ^ 0 as shown in 
Fig. 2(left). This leads to a rearrangement of the hole 
gas and the lowest energy subbands, and, in turn, the 
appearance of a kink in the emission line of the ground e- 
h state in the region around 5 ^ 10 T. As shown in Fig. 4, 
just at the complete depletion of the n = 3 hole level 
the gas redistribution affects the energy of the ground 
hole state at k = 0 which increases by ^4meV. Such an 
effect is related to the Fermi level pinning, that is 
goes from being pinned at the n = 3 subband to the 
flat n = 1,2 ones as found in other works.Note that 
in the present calculation PL spectra originate from the 
recombination of the hole and electron states with k = 0, 
since this is the only occupied state in the conduction 
band. At larger fields, when the dispersion flattens, also 
finite wavevectors are occupied. However, such states are 
localized on opposite sides of the structures for electrons 
and holes with the same k. Hence, the overlap decreases 
and, as a consequence, the PL intensity decreases as well. 



FIG. 3. Square modulus of the lowest /c = 0 electron (top) 
and hole (bottom) envelope function at selected fields, as in¬ 
dicated, for the sample of Fig. 2. 



B(T) 


FIG. 4. Energy and subband occupation for the sample of 
Fig. 2. Top: Energy of the ground hole state at k = 0 vs B. 
The horizontal line indicates the hole Eermi level. Bottom: 
Occupation of n = 1 (blue line) and n = 3 (cyan line) hole 
states at /c = 0 as a function of the magnetic field. 


In the high field regime, the emission energy of the 
e-h ground state increases linearly with the field, as ex¬ 
pected. The peak also reduces its intensity as a conse¬ 
quence of both the decrease of the overlap integral be¬ 
tween hole and electron states (see Fig. 3) and the de¬ 
creasing occupation of the electron ground state at k = 0. 

At = 10^ cm“^ the photo-excited hole density is of 
the same order of the free-charge density (8.3 x 10^ cm“^ 
at H = 0 T, 8.4 x 10^ cm“^ at H = 22 T). Figure 5 dis¬ 
plays the subband dispersion and the PL spectra for this 
configuration. The overall dispersion is similar to the 
previous case even though electrons now occupy a few 
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excited subbands at the lowest field. The PL spectrum 
at B = 0 exhibits four peaks, with the one due to e-h 
ground state at about 1520 meV being very weak. Addi¬ 
tional stronger peaks at 1523 meV, 1526 and 1530 meV 
are due to transitions between excited states. All emis¬ 
sion lines exhibit a diamagnetic shift at low fields and a 
linear one at high magnetic fields. At increasing fields, 
the number of emission lines decreases since the occupa¬ 
tion of the upper electron levels vanishes. At 5 = 10 
only two peaks are visible, one related to the e-h ground 


state and the other one originating from the recombina¬ 
tion of an electron in the ground state and a hole in the 
excited state (n = 3) spreading over the GaAs core (see 
the inset of the right panel of Fig. 5). As B increases, 
the n = 3 subband gets closer to the lowest one in the 
/c-interval around k = 0. This is due to the fact that 
the n = 3 hole state develops lobes localized on oppo¬ 
site sides of the structure where potential energy is lower 
with respect to the GaAs core. Therefore the two peaks 
approach each other and merge in an unique peak whose 
emission energy increases linearly with B. 
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FIG. 5. Same as Fig. 2, but with = 10® cm ^ 
n = 1 electron (left) and n = 1 (top right) and n 
B = 10 T. 


. The inset in the right panel displays the squared envelope functions for the 
= 3 hole (bottom right) states whose recombinations yield the two peaks at 


In the large excitation power regime, namely = 
10^ cm“^, the density of photoexcited holes is about 
equal to free hole density (0.85 x 10^ cm“^ at 5 = 0, 
0.84 X 10^ cm“^ at 5 = 22 T). Therefore, the electron 
density is now about one half the total density of holes, 
and a large number of both conduction and valence sub¬ 
bands are occupied, as shown in Fig. 6. Consequently, 
the low field PL spectrum shown in Fig. 6 is consider¬ 
ably more complex than in the previous cases, cover¬ 
ing the range 1523 meV to 1555 meV. As B rises, the 
number of peaks decreases and in the high-held region 
three PL emission lines can be identihed, each originat¬ 
ing from diagonal transitions involving electron and hole 
states with n = 1, 2, 3 (the square modulus of these wave- 
functions is reported in Fig. 7). Unlike the low-excitation 
power regime, the ground hole state also spreads over the 


central region due to the Coulomb attraction from free 
electrons. The peak at 1545 meV for B= 10 T related to 
the transition q ^ o decreases its intensity with the 
held while the other two emission lines due to the n = 1, 2 
diagonal transitions increase unlike what observed before 
in the PL spectra for lower excitation power. Here, the 
occupation number of n = 1,2 electron and hole levels 
at /c = 0 is always equal to 1 for any value of H, and 
the intensity of the emission lines is only modulated by 
the overlap integral which increases at larger helds. The 
lowest electron and hole states exhibit the same localiza¬ 
tion: the ground states spread over the core while the 
n = 2 ones are concentrated in two lobes localized on 
opposite sides. Such features are more marked for larger 
helds and, thus the integral overlap between the electron 
and holes states rises. 
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FIG. 6. Same as Fig. 2, but with = 10^ cm 



FIG. 7. Square modulus of the lowest k = 0 electron (top) 
and hole (bottom) envelope function at S = 22 T for the 
sample of Fig. 6. 


2. n-doping 

We next study the fan diagrams obtained in a n- 
doped sample with no = 4.5 x 10^^ cm“^ and a = 
10^ cm“^ with the same field configuration as in the 
previous section. Now the hole density is only due to 
photo-excitation, while the electron density due to photo¬ 
excitation is small, but not negligible, with respect to the 
free charge induced by n-doping (9.4x 10^ cm“^ at 5 = 0, 
9.5 X 10^ cm“^ at 5 = 22 T). From this point of view, 
this regime corresponds to the second case analyzed in 
p-doped systems (and reported in Fig. 5). Due to the 
large donor density, which induces an inversion layer at 


the hetero-interface, in the n-doped case the electron sub¬ 
bands deviate significantly from parabolicity, while holes 
tend to form Landau-like and edge states. For sufficiently 
large fields, the electron subbands develop minima which 
correspond to states localized at the flanks. 

The right panel of Fig. 8 displays the corresponding 
PL spectra. For B = 0 we observe several emission 
lines. Note that the lowest one is at 1515 meV, below 
the bulk GaAs band gap. Indeed, as shown in Fig. 9 (top 
panel), the hole ground state spreads over the GaAs core 
while the electron ground state is localized at the inner 
GaAs/AlAs interfaces depleting the central zone. The 
bending of the bands due to electrostatic potential mostly 
generated by the dopants is more pronounced in the re- 



























FIG. 8. Same as in Fig. 2 but for a n-doped sample with with ud = 4.5 x 10 ^ and nf 
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gions farther away from the central core. Therefore, the 
down shift of electron levels is generally larger than for 
hole states, and, as a consequence, the energy of ground 
e-h state lies below the bulk GaAs band gap. Such a be¬ 
havior was not observed in the p-doping regime since in 
that case the band bending was less pronounced due to 
the smaller density of dopants {ua = 2.5 x 10^^ cm“^ vs 
no = 4.5 X 10^^ cm“^). 

In this sample, the number of transitions at low fields 
is larger than a for p-doped sample in the corresponding 
regime of photo-excitation power, although the number 
of subbands occupied by the minority carriers is simi¬ 
lar. This is due to the different localization of the ma¬ 
jority charge which, for n-type doping (although mainly 
localized at the hetero-interface, as seen in the bottom 
panel of Fig. 9), extends to the core and overlaps with 
the hole ground state. As the magnetic field increases, 
the number of emission lines decreases due to hole state 
depopulation. At 5 = 10 T, only the hole ground state 
and few electron subbands are occupied while the Fermi 
level is pinned, as shown in left panels of Fig. 8. In the 
intermediate fields regime we identify three main peaks 
originating from the recombination of the hole ground 
state with different electron levels. 

The emission energy due to the e-h ground state ex¬ 


The spectra exhibit the same qualitative behavior, with 
several emission lines at low B which reduce to one peak 
(ground e-h state) at larger fields. However, the emission 
intensity is anisotropic, especially at high field intensity. 
At H = 22 T, the peak intensity at 0 = tt/S is about 
twice as strong as the one at 0 = 7r/2. This is due to 
the different localization of the n = 1 electron and hole 
states at /c = 0, shown in the right panel of Fig. 11. The 
hole state shows two lobes localized by the strong carrier- 


hibits a fiat region between B = 9 and 12 T. For B > 9T 
the occupation of the n = 1 hole state at k = 9 decreases 
significantly since the lowest subband flattens and states 
with k 0 are increasingly populated at the expense of 
the k = 0 state, as shown in the top panels of Fig. 10. 
Since finite k states are localized away from the NW axis, 
the field induces a restructuring of the hole density which 
becomes anisotropic (see the bottom panel of Fig. 10). 

For the high-held regime, only one peak survives, due 
to the recombination of the ground e-h state. 

C. Anisotropy 

Carriers are either homogeneously distributed around 
the hetero-interface or form ID and 2D channels accord¬ 
ing to the device parameters. In the latter case, we expect 
the direction of the held to have a remarkable impact. 
Therefore, we next compare PL spectra with the held 
oriented in two non-equivalent directions, namely per¬ 
pendicular to the top/bottom facets {0 = 7r/2) and along 
a maximal diameter {0 = tt/S). We assume p-doping 
with Ua = 2.5 x 10^^ cm“^ at the intermediate excita¬ 
tion power regime, that is = 10^ cm“^. Left panels 
of Fig. 11 display the PL spectra evaluated for the two 
conhgurations at diherent values of the held intensity. 


carrier repulsion towards the bottom and top facets (at 
0 = 7r/2) or corners (at 0 = tt/S). Since the facet-to- 
facet distance is the shorter, due to Coulomb interaction 
the energy of the state in the hrst conhguration turns 
out to be a few tenths of meV greater than the one with 
the held along the maximal diameter. As a consequence, 
its occupation is smaller. The electron states elongate in 
the held direction and, also in this case, the energy level 
obtained at 0 = 7 t/2 turns out to be greater than the one 
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n-doping p-doping 



FIG. 9. Top: self consistent conduction (black line) and va¬ 
lence (red line) band profile and the square modulus of the 
n = 1 electron (green dashed line) and hole (blue dotted line) 
wavefunction along the y-axis at B = 0. Bottom: electron and 
hole charge densities with — 10® cm“^ at B = 2 T for 
riD = 4.5 X 10^^ cm“^ (left panels) and ua = 2.5 x 10^^ cm“^ 
(right panels). 


found at 0 = tt/S, and this yields a discrepancy in the 
occupation of the states. Thus, the higher population of 
both n = 1 hole and electron states leads to the stronger 
intensity of the ground e-h state emission line at 0 = tt/3 
with respect to the one found at 0 = 7r/2. 


IV. SUMMARY AND CONCLUSIONS 

We have investigated theoretically the magneto-optical 
properties of hole and electron gases in GaAs/AlAs core- 
multishell NWs which are either n- or p-doped and sub¬ 
ject to a transverse magnetic field. The tailoring of PL 
spectra by the field has been rationalized in different dop¬ 
ing regimes. Calculations were performed by means of a 
self-consistent mean-field approach, with the inclusion of 
photogenerated hole-electron pairs that add to the free 




B = 9J B=11T B = 13T 



FIG. 10. Top left: n = 1 hole subband vs kl% evaluated for B 
ranging from 9 to 13 T. Top right: occupation of n = 1 hole 
state at k = 0 vs B. Bottom: hole gas density at selected 
fields, as indicated. 


carriers induced by the modulation doping. 

Our results indicate that the behavior of magneto-PL 
spectra offers a valuable guideline to estimate the density 
of the doping in CMS NWs. In fact, when the density 
of free carriers due to the doping is greater than the one 
of photogenerated charges, the spectra are quite simple 
and only the emission line due to the ground e-h state 
survives at high fields. Differently, when the excitation 
power is larger and the density of photoexcited carri¬ 
ers is comparable or higher with respect to the doping- 
induced free charge, several PL peaks are found, that 
can be ascribed to recombinations from specific excited 
states, whose character is, in turn, strongly modified by 
the magnetic field. 

We found a remarkable difference between samples 
with either electron or holes as majority carriers. Indeed, 
at low magnetic fields and low photoexcitation power, the 
qualitative behavior of PL spectra for n- and p-doped 
sample is quite similar although, in the case of n-doping, 
the number of emissions lines is greater than the one 
found for the p-doped sample. This is a consequence of 
the different localization pattern of the majority carri¬ 
ers. Still in the low-field regime, when the density of 
photo-created hole-electron pairs is significant, the emis¬ 
sion lines exhibit a diamagnetic shift, in both samples. 
This shift becomes linear for higher magnetic fields. How¬ 
ever, the main peak, due to the ground e-h state, devi¬ 
ates from the expected behavior at low values of density 
of photogenerated hole-electron pairs. Specifically, in the 
intermediate field regime, it exhibits a kink (p-doping) or 
a fiat region (n-doping). In both cases this is due to the 
rearrangement of the gas of holes, whose heavier effec¬ 
tive mass eases their localization near the interface with 
respect to the electrons. In the high power excitation 
regime, the large Coulomb interaction prevents the hole 
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B = 2 T B = 10 T B =22 T 




FIG. 11. Left: PL spectra in different directions (top) and absolute difference between the PL spectra (bottom) at selected 
fields, as indicated. Right: square modulus of the n = 1 hole and electron envelope functions at /c = 0 and R = 22 T for the 
two field orientations. 


gas rearrangement. 

The differences in the magneto-PL spectra illustrated 
in this work are a valuable tool to discriminate not only 
the free-carrier density regime and type, but also the 
relative density of doping- and photoexcitation-induced 
carriers. This, together with the assessment of the free- 
carrier spatial distribution, is pivotal in the characteriza¬ 
tion and engineering of CMS NWs. 


ACKNOWLEDGMENTS 

Numerical simulations were performed at CINECA 
within the Iscra C project MPL-CSNW. We acknowledge 
partial financial support from University of Modena and 
Reggio Emilia, with Grant “Nano- and emerging materi¬ 
als and systems for sustainable technologies.” M. R. ac¬ 
knowledges financial support from APOSTD/2013/052 
Generalitat Valenciana Vali+d Grant,. 


* fabrizio.buscemi@unimore.it 

^ G. M. Lieber and Z. L. Wang, MRS Bulletin 32 , 99 (2007). 

2 Y. Gui, Z. Zhong, D. Wang, W. U. Wang, and G. M. 
Lieber, Nano Letters 3 , 149 (2003). 

^ J. A. Gzaban, D. A. Thompson, and R. R. LaPierre, Nano 
Letters 9, 148 (2009). 

^ M. J. Holmes, K. Ghoi, S. Kako, M. Arita, and 
Y. Arakawa, Nano Letters 14, 982 (2014). 

^ F. Qian, S. Gradeak, Y. Li, G.-Y. Wen, and G. M. Lieber, 
Nano Letters 5, 2287 (2005). 

® H. Shtrikman, R. Popovitz-Biro, A. Kretinin, and 
M. Heiblum, Nano Letters 9, 215 (2009). 

^ H. Shtrikman, R. Popovitz-Biro, A. Kretinin, L. Houben, 
M. Heiblum, M. Bukaa, M. Galicka, R. Buczko, and 
P. Kacman, Nano Letters 9, 1506 (2009). 

® K. Tomioka, J. Motohisa, S. Hara, K. Hiruma, and 
T. Fukui, Nano Letters 10 , 1639 (2010). 

^ M. J. Tambe, S. K. Lim, M. J. Smith, L. F. Allard, and 
S. Gradeak, Applied Physics Letters 93, 151917 (2008). 

S. Funk, M. Royo, 1. Zardo, D. Rudolph, S. Morktter, 
B. Mayer, J. Becker, A. Bechtold, S. Matich, M. Dblinger, 
M. Bichler, G. Koblmller, J. J. Finley, A. Bertoni, 


G. Goldoni, and G. Abstreiter, Nano Letters 13 , 6189 
(2013). 

M. Fickenscher, T. Shi, H. E. Jackson, L. M. Smith, J. M. 
Yarrison-Rice, G. Zheng, P. Miller, J. Etheridge, B. M. 
Wong, Q. Gao, S. Deshpande, H. H. Tan, and G. Jagadish, 
Nano Letters 13 , 1016 (2013). 

J. Jadczak, P. Plochocka, A. Mitioglu, 1. Breslavetz, 
M. Royo, A. Bertoni, G. Goldoni, T. Smolenski, P. Kos- 
sacki, A. Kretinin, H. Shtrikman, and D. K. Maude, Nano 
Letters 14, 2807 (2014). 

T. Shi, H. E. Jackson, L. M. Smith, N. Jiang, Q. Gao, 

H. H. Tan, G. Jagadish, G. Zheng, and J. Etheridge, Nano 
Letters 15 , 1876 (2015). 

B. M. Wong, F. Lonard, Q. Li, and G. T. Wang, Nano 
Letters 11 , 3074 (2011). 

A. Bertoni, M. Royo, F. Mahawish, and G. Goldoni, Phys. 
Rev. B 84, 205323 (2011). 

K. Tomioka, M. Yoshimura, and T. Fukui, Nature 488, 
189 (2012). 

H. Ajiki and T. Ando, Journal of the Physical Society of 
Japan 62 , 1255 (1993). 

G. Ferrari, A. Bertoni, G. Goldoni, and E. Molinari, Phys. 




































11 


Rev. B 78 , 115326 (2008). 

S. Bellucci and P. Onorato, Phys. Rev. B 82 , 205305 

( 2010 ). 

G. Ferrari, G. Goldoni, A. Bertoni, G. Guoghi, and 
E. Molinari, Nano Letters 9 , 1631 (2009). 

M. Royo, A. Bertoni, and G. Goldoni, Phys. Rev. B 87 , 
115316 (2013). 

M. Royo, G. Segarra, A. Bertoni, G. Goldoni, and 
J. Planelles, Phys. Rev. B 91 , 115440 (2015). 

L. V. Titova, T. B. Hoang, J. M. Yarrison-Rice, H. E. Jack- 
son, Y. Kim, H. J. Joyce, Q. Gao, H. H. Tan, G. Jagadish, 
X. Zhang, J. Zou, and L. M. Smith, Nano Letters 7 , 3383 
(2007). 

T. B. Hoang, L. V. Titova, J. M. Yarrison-Rice, H. E. 
Jackson, A. O. Govorov, Y. Kim, H. J. Joyce, H. H. Tan, 
G. Jagadish, and L. M. Smith, Nano Letters 7 , 588 (2007). 
J. L. Boland, S. Gonesa-Boj, P. Parkinson, G. Ttncoglu, 
E. Matteini, D. Rffer, A. Gasadei, E. Amaduzzi, E. Jabeen, 
G. L. Davies, H. J. Joyce, L. M. Herz, A. Eontcuberta i 
Morral, and M. B. Johnston, Nano Letters 15 , 1336 
(2015). 

L. M. Smith, H. E. Jackson, J. M. Yarrison-Rice, and 
G. Jagadish, Semiconductor Science and Technology 25 , 
024010 (2010). 

P. Plochocka, A. A. Mitioglu, D. K. Maude, G. L. J. A. 
Rikken, A. Granados del guila, P. G. M. Ghristianen, 


P. Kacman, and H. Shtrikman, Nano Letters 13 , 2442 
(2013). 

The small Zeeman energy (<0.3 meV at the largest field 
used) is below the thermal linewidth^^ and has been ne¬ 
glected. 

E. Rossi and E. Molinari, Phys. Rev. Lett. 76 , 3642 (1996). 
S. Sedlmaier, M. Stopa, G. Schedelbeck, W. Wegscheider, 
and G. Abstreiter, Phys. Rev. B 65 , 201304 (2002). 

We do not include the exchange and correlations effects in 
the mean field since they turn out to be negligible. 

Due to the symmetry of the Hamiltonian in Eqs. (2), states 
with negative k can be obtained from (r, k) (—r, —k). 
G. Bastard, Wave mechanics applied to semiconductor het¬ 
erostructures, Monographies de physique (Les Editions de 
Physique, 1988). 

M. D. Pashley, K. W. Haberern, R. M. Eeenstra, and P. D. 
Kirchner, Phys. Rev. B 48, 4612 (1993). 

M. Levinshtein, S. Rumyantsev, and M. Shur, Hand¬ 
book Series on Semiconductor Parameters (World Scien¬ 
tific Publishing, 1996). 

G. Eishman, Phys. Rev. B 52 , 11132 (1995). 

G. Goldoni, E. Rossi, E. Molinari, and A. Fasolino, Phys. 
Rev. B 55, 7110 (1997). 

P. YU and M. Gardona, Fundamentals of Semiconductors: 
Physics and Materials Properties (Springer, 2013). 

B. Jogai, Journal of Applied Physics 91 , 3721 (2002). 



